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ABSTRACT 

A candidate diffuse stellar substructure was previously reported in the halo of the nearby dwarf star- 
burst galaxy NGC 4449 by Karachentsev et al. We map and analyze this feature using a unique com- 
bination of deep integrated-light images from the Black Bird 0.5-meter telescope, and high-resolution 
wide-field images from the 8-meter Subaru telescope, which resolve the nebulosity into a stream of red 
giant branch stars, and confirm its physical association with NGC 4449. The properties of the stream 
imply a massive dwarf spheroidal progenitor, which after complete disruption will deposit an amount 
of stellar mass that is comparable to the existing stellar halo of the main galaxy. The ratio between 
luminosity or stellar- mass between the two galaxies is ^ 1:50, while the indirectly measured dynam- 
ical mass-ratio, when including dark matter, may be ~ 1:10-1:5. This system may thus represent a 
"stealth" merger, where an infalling satellite galaxy is nearly undetectable by conventional means, yet 
has a substantial dynamical influence on its host galaxy. This singular discovery also suggests that 
satellite accretion can play a significant role in building up the stellar halos of low-mass galaxies, and 
possibly in triggering their starbursts. 



L INTRODUCTION 

A fundamental characteristi c of the modern cold dark 
matter (ACDM) cosmology (|Mo et all 120101 ) is that 
galaxies assemble hierarchically under the influence of 
gravity, continually accreting smaller DM halos up un- 
til the present day. If those halos contain stars, then 
they will be visible as satellites around their host galax- 
ies, eventually disrupting through tidal forces into dis- 
tinct streams_a£d shells before phase- mixing into obscu- 
rity (jCooper et al.ll20l"ol ). This picture seems to explain 
the existence of satellites and substructures observed 
around massive galaxies ( Arp 1966; Schwcizcr & Scitzer 
1988t iMcConnachie et all 120091: iMartinez-Delgado et all 
2010f ). However, quantitative confirmation of this as- 
pect of ACDM has been more elusive, with linger- 
ing doubts provoked by small-scale substructure obser- 
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vations (iL ovell et all [20121 : iBovlan-Kolchin et all [20121 : 
iFerrero et al. 2012). 

There has been relatively little work on substruc- 
ture and merging in the halos around low-mass, 
"dwarf" galaxies. Many cases of extended stellar ha- 
los around dwarfs have been identified observationally 
(jStinson et al.1120091 ). but it is not clear if these stars were 
accreted, or formed in-situ. Star formation in d warfs is 
thought to occur in st ochastic episodes (Tolsto y et al.l 
l2009tlWeisz et al.ll2011[ ). which could be triggered by ac- 
cretion events. 

An iconic galaxy in this context is NGC 4449, a dwarf 
irregular in the field that has been studied intensively 
as one of the strongest galaxy-wide starbursts in the 
nearby universe. Its absolute magnitude of My = —18.6 
makes it an LMC-analogue (and not formally a dwarf by 
some definitions), but with a much higher star formation 
rate. It is strongly suspected to have recently interacted 
with another galaxy based on various signatures includ- 
ing peculiar kinematics in i ts cold g as and HII regions 
(iHartmann et al. 1986; Hun ter et al.i ri998). but the na- 
ture of this interaction is unknown. 

An elongated dwarf galaxy or str eam candidate near 
NGC 4449 was first noticed by iKarachentsev et"an 
(pUol from Digitized Sky Survey (POSS-II) plates (ob- 
ject dl228+4358), and is visible in the Sloan Digital Sky 
Survey (SDSSJI3. Here we present deep, wide-field opti- 
cal imaging that supplies the definitive detection of this 
ongoing accretion event involving a smaller galaxy, lead- 
ing to interesting implications about the evolution of this 
system and of dwarf galaxies in general. 



http:/ /www. sdss.org/ 
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2. OBSERVATIONS AND DATA REDUCTION 

Our observations of NGC 4449 and its surroundings 
consist of two main components. The first is imaging 
from a small robotic telescope, to confirm the presence 
of a low-surface-brightness substructure and provide its 
basic characteristi cs (similar techniques were us e d wit h 
larger galaxies in iMartmez-Delgado et alJ I2008L l2010t ). 
The second is follow-up imaging with the Subaru tele- 
scope to map out the resolved stellar populations. 

We obtained very deep images with the f/8.3 Ritchey- 
Chretien 0.5-meter telescope of the Black Bird Remote 
Observatory fBBROf^ during different dark-sky observ- 
ing runs over the periods 2010-04-13 through 2010-06-10, 
and 2011-01-13 through 2011-01-28 (UT). We used a 16 
mega-pixel Apogee Imaging Systems U16M CCD cam- 
era, with 31.3' X 31.3' field-of-view and 0.46 arcsec pix~^ 
plate-scale. We acquired 18 hours of imaging data in 
half-hour sub-exposures, using a non-infrared clear lumi- 
nance (A = 3500-8500A) Astrodon E-series filter. Each 
sub-exposure was reduced following standard procedures 
for dark-subtraction, bias-cor rection, and flat-fielding 
(|Martmez-Delgado et aT]|2009D . 

The resulting image was calibrated photometrically 
to SPSS using the brighter regions of NGC 4449 (see 
IMartmez-Delgado et al.ir2010t) . The final image has 5-a 
(jf-band surface-brightness detection limits from 26.4 to 
27.5 mag arcsec"^ for seeing-limited and large-scale dif- 
fuse features, respectively. 

We subsequently obtained images from the 8.2-m Sub- 
aru Telescope and the Suprime-Cam wide-field imager 
(34'x 27' field-of-view, 0.202" pixel-scale; Mivazaki et al, 
I2002D on 2011-01-05 (UT). Conditions were photometric, 
and we took dithered exposures in r' and i' bands, with 
total exposure times of 225 s per filter. We reduced the 
data using a modified SDFRED pipeline (jOuchi et al.l 
120041 ) ■ including bias subtraction, fiat-fielding, and dis- 
tortion correction. Each frame was re-projected to a 
common astrometric coordinate system followed by back- 
ground rectification and image co-addition using Mon- 
tageES. 

The exquisite image quality 0.5" FWHM) al- 
lows us to resolve individual stars in the outer regions 
of NGC 4449. We carried out point-spread- function 
photo metry using DAOPHOT II/ALLSTARS (jStetsonl 
Il987f) . and identified stars as objects with sharpness 
parameter \S\ < 1.0. We calibrated the photometry 
based on two central images from the Hubble Space 
Telescope Advanced Camera for Surveys (HST/ACS) 
(jAnnibali et al .' 2008, hereafter A+08). 

The ACS photometry was originally in F555W 
and F814W, and recalibrated to Johnson-Cousins VI 
(A-l-08). We used fairly bright, red stars in common be- 
tween the datasets to derive linear transformation equa- 
tions from r'i' instrumental magnitudes to VI, including 
foreground extinctio n corrections of E{B — V) — 0.019 
(jSchlegel et al.|[r998D . Our final star catalog has statisti- 
cal internal errors inV — I color of ~0.11 mag and ~0.14 
mag at V=25 and V=2Q, respectively. 



BBRO was originally situated in the Sacramento Mountains 
(New Mexico, USA), and later moved to the Sierra Nevada Moun- 
tains (California, USA). 

http:/ /montage. ipac.caltech.edu/ 



3. STREAM MORPHOLOGY 

Figure [l](te/it) shows a BBRO image subsection, where 
with an adopted distance of 3.82 Mpc (A4-08), 1' corre- 
sponds to 1.1 kpc. Clearly visible near the minor axis 
of NGC 4449, ~ 10 kpc to the southeast, is a very elon- 
gated, S-shaped feature of dimensions ~ 1.5 x 7 kpc, 
which we designate "the stream" . The stream's position 
does not overlap with any o f the compl ex Hl-gas fea- 
tures surrounding NGC 4449 (|Hunter et al . 1998). Also, 
it is on the opposite side of the galaxy with respect to 
an interesting star clu ster that may be linke d to another 
past accretion event (jAnnibali et al.|[2012|) . The main 
galaxy's existing stellar halo is also apparent, including 
shell-like features in the southwest previously noticed by 
iHunter et aP (|1999f ). 

The right-hand panel shows a Suprime-Cam zoom-in 
of the stream. It is clearly resolved into stars and has the 
general appearance of a dwarf spheroidal (dSph) galaxy 
that is elongated by tides. This classification (defined 
by a spheroidal morphology, stellar mass < 10^ Mq, 
and the absence of star formation an d cold gas) is also 
suppo rted by the HI non-detection of iHuchtmeier et al.l 

mm . 

To complement these integrated surface-brightness 
maps, we construct stellar-density maps around 
NGC 4449 using individual RGB stars down to g' = 25.8 
from Suprime-Cam. We discuss the RGB selection and 
modeling later, but in general these stars trace a popula- 
tion older than 1 Gyr. We define a grid across the image 
with 130 X 130 pc^ bins, and count the number of stars 
in each bin, subsequently applying a Gaussian smooth- 
ing kernel with (Tsmooth = 130 pc. There are typically 
7-9 stars per bin in the stream. 

Figure [2][a) shows the resulting Subaru-based stellar- 
density map, where the overall stream morphology is in- 
distinguishable from the BBRO integrated-light results. 
Since there is no evidence for resolved young stars (see 
next section), we infer that the visible light is dominated 
by RGB stars and hence by a population older than 1 
Gyr. 

Panel (b) shows a higher-contrast version of the same 
map, demonstrating that the main galaxy's pre-existing 
stellar halo extends out to at least '^10 kpc. We can 
furthermore discern a very faint feature that seems to 
extend the stream's angular path in a loop-like struc- 
ture. This loop is also apparent in the BBRO image, 
and we infer that the stream is stretched out over at 
least half its orbit, with the projected turning point at a 
galactocentric radius of 13 kpc. 

We next wish to locate the disrupting satellite galaxy's 
original center or nucleus. In both the BBRO and Sub- 
aru images there is a density enhancement near the mid- 
point of the "S" , as would be expected if the "arms" are 
leading and trailing tidal tails around a still marginally- 
bound, or just disrupting, main body. We construct 
a stellar-density contour map for the stream region, 
using RGB stars and tJsmooth = 270 pc. The cen- 
tral stellar clump is visible in Figure [2][d) , with a po- 
sition slightly offset from the stream's main ridgeline 
(aj20oo = 12 28 43.32, Jjsooo = +43 58 30.0; positional 
uncertainty ~ 6"). An off-center nucleus is also seen 
in a tidally-disrupting Milky Way dSph, Ursa Major 
(|Martfnez-Delgado et aT]|2001[ ). 



Fig. 1.— NGC 4449 and its halo stream. Left: image from BBRO, showing a 19.0' X 24.5' (21 X 27 kpc) field. Right: 5.5' X 8.6' 
(6 X 9.5 kpc) subsection of the Subaru/Suprime-Cam data, showing the stream resolved into stars. In both panels, shallower BBRO 
exposures in red/green/blue filters provide indicative colors. 



4. STELLAR POPULATIONS 

Figure [3] shows the color-magnitude diagram (CMD) 
for point sources in the stream region. The RGB 
stars are the dominant feature, along with a few 
brighter, redder stars that may be oxygen- or carbon-rich 
thermally-pulsating asymptotic giant bran ch stars from 
an in termediate-age or old population (iMarigo et al.l 
120031 ). We find no blue stars that would trace recent 
star formation. 

The detection of the tip of the RGB (TRGB) per- 
mits a_distaj]£e estimate. Using techniques from A-f08 
and iCioni etaHl (poM l. we find a TRGB magnitude of 
^TRGB = 24.06 ± 0.04 (random) ±0.08 (systematic). 
The random error was estimated using bootstrapping 
techniques; the systematic error is dominated by the 
magnitude-transformation uncertainties. The main body 
of NGC 4449 was found by A-l-08 to have /trgb = 
24.00 ± 0.01 (random) ±0.04 (systematic). Thus the 
stream is at the same distance as the main body, to 
within 180 kpc, and we conclude that there is a phys- 
ical association rather than a chance superposition. 

Although the RGB is affected by the well-known age- 
metallicity degeneracy, this feature can still be used to 
constrain the properties of stars older tha n ^ 1 Gyr. In 
FigureOwe overplot the Padua isochrones fG irardi et al.l 
[2002 ) for ages of 2, 4, and 10 Gyr, for both Z = 0.004 
and Z = 0.001. For Z = 0.004 the 2-4 Gyr isochrones 



trace the data reasonably well. The Z = 0.001 models 
appear significantly bluer than the mean RGB color (by 
> 0.13 mag: our systematic color uncertainties are ^0.1 
mag), although a Z — 0.001, 10 Gyr model is consistent 
with the blue edge of the observed RGB. Presumably a 
10 Gyr model with Z 0.002 would also be consistent 
with the data, while perhaps providing a better match to 
the observed CMD slope. We conclude that, if the bulk 
of the RGB stars are old (> 10 Gyr), the metallicity 
range is roughly Z = 0.001-0.004, while for younger ages 
the metallicity range shifts to higher values. 

These results are comparable to the RGB analysis 
of the main body of NGC 4 449 by A+08 (their fig- 
ure 17; see also lRvs et al.l[20T]l ). Therefore both the main 
galaxy and the stream contain similar old, intermediate- 
metallicity populations, although the main galaxy also 
contains very young stars, as well as mo re metal-poor 
stars as inferred from its globular clusters (jStrader et al.l 
I20TI . 

We provide a preliminary overview of spatial stellar 
population variations by splitting the RGBs into color- 
based subpopulations, using the 4 Gyr, Z — 0.004 model 
as a boundary. We then create stellar-density maps as 
before, for the two subpopulations separately. Using blue 
and green color-coding to represent the subpopulations 
left and right of the model boundary, we show the results 
in Figure [2l^c). The stream's bright parts have no visible 
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Fig. 2. — Stellar-density maps of NGC4449 and its stream, based on RGB counts, (a) 30' X 30' (33 X 33 kpc) image with linear scaling. 
A compact density enhancement within the stream may be the progenitor galaxy's remnant nucleus, (b) The same data with less dynamic 
range. The stream shows a loop-like structure curving back toward the main galaxy, which in turn also shows a shell-like overdensity in 
its halo toward the southwest, reaching similar projected galactocentric radii as the stream, (c) Composite map of the stream, color-coded 
blue for bluer RGB colors, and green for redder colors (see main text). The RGBs in the stream's main part have similar average colors to 
the main galaxy's inner halo stars, with redder populations apparent on the outskirts of both the halo and stream, (d) Star-count contour 
map of the stream, with contour levels corresponding to 5-10 stars per 50 arcsec^ bin, in intervals of 1 star per bin. 



RGB color gradient, and have an overall color similar to 
the main galaxy's halo at radii of ~ 3-5 kpc. Both the 
stream's faint-loop continuation, and the halo at ~ 5- 
10 kpc, are redder, implying older or more metal-rich 
stars. 

5. STELLAR AND DYNAMICAL MASS 

We now proceed to estimate the luminosities and 
masses of NGC 4449 and its stream. For NGC 4449, 
it is straightforward to add up the SDSS pixel fluxes 
inside the "optical radius" {fir — 25 mag arcsec"^). 
Wc find an extinction-corrected Air = —17.8. For the 
stream, we use the SDSS-calibrated BBRO image, inte- 
grating the flux within the faintest isophote that closes 
without including the main galaxy, which corresponds to 
/ig = 26.75 mag arcsec"^ and a stream semi-major axis 
distance of 3 kpc. We find a stream magnitude of Mr = 
— 13.5, which is comparable to the brightest Local Group 
dSphs, Fornax and And VII. Such galaxies have typical 
proje cted half-light radii of ^ 0.4-1.0 kpc (jBrodie et al.l 
1201 If ), which is consistent with the stream's ^ 0.8 kpc 
half-width. 

For both galaxies, these luminosity estimates are lower- 
limits since they do not include potential extended low- 
surface-brightness features. The implied luminosity-ratio 
is - 1:50. 

To calculate stellar ma sses M-^, we us e two differ- 
ent approaches, adopting a lChabrieil ()2003f ) initial-mass- 
function (IMF; final mass, including stellar remnants). 
The first is based on the integrated optical colors, follow- 
ing the relations between color and stellar mass-to-light- 
ratio (T,^) from iZibetti et all ([2009L Table Bl). This 
paper also introduced a technique for mapping out lo- 
cal Ti, and M^, variations pixel-by-pixel, which we ap- 
ply to NGC 4449, and after integrating, find a total 
Mi, = 7.46 X IO^Mq. For the stream, we assume a 



global color corresponding to the central SDSS measure- 
ment, g — r — 0.45 ± 0.1. We find a stream mass of 
Mi, = l.S^Q g X 1O^M0, implying a stellar mass-ratio be- 
tween stream and host of ~1:50, for any uniform IMFFI 

The second mass-estimation approach uses the CMD, 
comparing observed star-counts to predicted numbers 
from a stellar populations model. For the stream, we use 
the /-band stellar luminosity function near the TRGB, 
and normalize it to Monte Carlo simulations drawn as- 
suming Z =0.001-0.004, and ages 2-10 Gyr. We obtain 
Mi, ^ (2-5.5) X 10'' Mq for the stream (a lower-limit be- 
cause of incompleteness), which agrees with the color- 
based results. For the main galaxy, a similar approach 
was followed bv lMcQuinn et all (|201Clf l , whose results im- 
ply Mi, = (1.2 ± 0.2) X IO^Mq. This mass is somewhat 
higher than by using colors, but the CMD-based stellar- 
mass-ratio comes out to be similar, ~ 1:40. 

Remarkably, the mass of NGC 4449 's stellar halo is 
si milar to th e strea m's mass: based on the RGB counts 
of iRvs et al.l (|2011[ ) and their normalization to K-hand 
surface-brightness photometry, we estimate Mi, ~ 2 x 
IO^A/q for the halo between projected radii of 5-10 kpc. 
This halo could have therefore been built up directly by 
one or a few past accretion events similar to the present- 
day stream. 

We next consider the dynamical masses of the host 
galaxy and its stream, including DM. The quantity that 
is arguably the most relevant to the current stream- 
galaxy interaction is the dynamical mass-ratio within 
the interaction region: the ~ 14 kpc galactocentric ra- 
dius. Based on the HI gas kinematics, we estimate an 
inclination-corrected circular velocity of Vc — 62 km s~^ 

The similarity of the luminosity- and mass-ratios implies that 
the luminosity-weighted estimate of T* for NGC 4449 happened 
to turn out the same as inferred for the stream. 
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Fig. 3. — Color-magnitude diagram of the stream region, centered 
(4.6 X 10.0 kpc) ficld-of-vicw, using r'i' photometry transformed to 
curves) and Z = 0.004 (red curves). For each Z, the isochrones show 

at this radius ()Baiaia et all 119941: iHunter et a"Lll2002[ ). 
which means a dynamical mass for the main galaxy of 
Mdy„(r < 15 kpc) - 1-1 X W^ Mq. The HI gas mass 
is ~ lO^M© (|Hunter et al.lll998f ). so this region is DM- 
dominated. Note that the Vc and Mi, values together 
suggest that NGC 4449 is intermedi ate in mass to the 
LMC and SMC (cf. iBesla et al.|[2010l) . 

For the stream mass, we have no direct measure- 
ments, and instead turn to a plausibihty argument 
based on Local Group dSphs, where the brightest 
cases have estimated circular ve locities of Vc 15- 
20 km s~^ on ~ 1-3 kpc scales ([Walker fc Penarrubial 
120111: iBovlan-Kolchin et all [20121) . Then if we assume 
the Vc values for both stream and main galaxy are fairly 
constant with radius, the ratio of yields the dynam- 
ical mass-ratio. This ratio is ^ 1:20-1:10, and thus the 
stream may be significantly perturbing the main galaxy. 

A final metric is the ratio of total (virial) halo masses 
Mvir, which are not directly measurable but may be in- 
ferred on a statistical basis, assuming a ACDM frame- 
work. In this context, it is well-established that the 
total mass-to-light-ratios of dwarf galaxies increase dra- 
matically at lower luminosities. Current estimates of 
Mi,-My\-c and luminosity-M yi]. relations (|Moster et al.l 
[2010t iToUerud et all [201ll) would imply Mvir ~ (1- 
5) X IO^Mq for NGC 4449, and a pre-infall mass of - (1- 
10) X 10"'^'' M0 for the stream progenitor — which although 
very uncertain, plausibly implies an initial virial mass- 
ratio of - 1:10-1:5. 

Thus what appears to be a very minor merger in visi- 
ble light may actually be closer to a major merger when 
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at (aj2000 = 12 28 43,(5j2ooo = +43 58 30.0), with a 250" x 540" 
VI. RGB isochrones are overlaid for metallicities Z = 0.001 (blue 
ages of 2, 4, and 10 Gyr, from left to right, respectively. 

including DM. Such an extreme circumstance could be 
compared with models of satellite disruption and poten- 
tiall y discriminate between A CDM and alternative theo- 
ries (IMcGaugh fc Woia[20T0l) . 

6. DISCUSSION 

We have detected and analyzed a stellar tidal stream 
in the halo of NGC 4449 which we interpret as the on- 
going disruption of a dSph galaxy by a larger dwarf (an 
LMC/SMC analogue^). This appears to be the lowest- 
mass primary galaxy with a verified stellar stream. 

We suggest some implications for galaxy evolution. It 
has been proposed that dSph's orbiting massive galaxies 
such as the Milky Way were "pre-processed" from gas- 
rich dwarfs by t i dal e ffects within dwarf-galaxy groups 
(ID'Onghia et al. 1 120091) . We may be witnessing such a 
transformation in-action, with the HI streams surround- 
ing NGC 4449 representing additional tidal debris. 

We also suspect it is not just a coincidence that such 
a novel stream was found first around one of the most 
intensely star-forming nearby galaxies. The accretion 
event may well be the starburst trigger. The period of el- 
evat ed star formation app ears to have started ~ 0.5 Gyr 
ago (jMcQuinn et al.l[20100 . which is suggestively similar 
to the stream's ~ 1-2 Gyr orbital perioc0 (and to any 
process that is linked to the dynamical time on ^ 30 kpc 

The LMC and SMC may also have a history o f interaction, 
with a stellar mass-ratio of ~ 1:15 IBesla et al.ll20l"0t) . 

Given a projected apocentric radius of i?a = 13 kpc, 
a circular orbit provides a lower-limit for the period of T = 
27riJa/nc ~ 1.3 Gyr. 
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scales). 

Are such accretion events frequent among other dwarf 
galaxies in recent epochs? We suspect that exact ana- 
logues to this stream are not very common, or they would 
have been noticed already in DSS/SDSS images. How- 
ever, if the stream had been only a bit fainter, more 
diffuse, or at a larger radius, it could have been missed, 
and thus there may be many more dwarf-hosted stellar 
streams awaiting detection. 

In theory, the history of DM halo assembly should be 
fairly scale- free, and ~ 1:10 mergers are expected to be 
the most generally d ominant contributors to mass growth 
(jStewart et al.ll2008i) . It is also increasingly recognized 
that such relatively minor mergers can have important 
effects on the larger galaxies, such as inciting global disk 
instabilities (Purcell et al. 2011). 

If streams as in NGC 4449 are common in dwarfs, they 
re-ignite classic ideas about galaxy interactions trigger- 
ing starbursts. Given the high rates of star formation in 
dwarf galaxies, it is natural to ask if satellites are respon- 
sible. Surveys along these lines have pr oduced mixed 
resul ts (jNoeske et al.ll2001t iBrosch et all 120041: iLi et al.l 
120081 ). but until now, low-surfacc-brightness objects such 
as dSphs would have been missed. 



Regardless of the implications for starbursts, dSph ac- 
cretion appears to be an increasingly viable avenue for 
direct assembly of dwarf galaxies' stellar halos — as wit- 
nessed by NGC 4449, and by Fornax, whi ch shows traces 
of sw allowing an even smaller dSph (iColeman et al.l 
l2005f l. Future observational determinations of dwarf 
stream frequency in combination with theoretical models 
may provide clues to the general substructure problem. 
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